International Journal of Modern Physics D 
© World Scientific Publishing Company 



SPECIAL PURPOSE PULSAR TELESCOPE FOR THE 
DETECTION OF COSMIC GRAVITATIONAL WAVES 



Shou-Guan Wang 
National Astronomical Observatories, Chinese Academy of Sciences, 
Beijing 100012, China 

Zong-Hong Zhu* 

National Astronomical Observatories, Chinese Academy of Sciences, 
Beijing 100012, China 
National Astronomical Observatory, 
2-21-1, Osawa, Mitaka, Tokyo 181-8588, Japan 

Zhen-Long Zou 

National Astronomical Observatories, Chinese Academy of Sciences, 
Beijing 100012, China 

Yuan-Zhong Zhang 
Institute of Theoretical Physics, Chinese Academy of Sciences, 
Beijing 100080, China 



Received (received date) 
Revised (revised date) 

Pulsars can be used to search for stochastic backgrounds of gravitational waves of cosmo- 
logical origin within the very low frequency band (VLF), 10~ 7 to 10~ 9 Hz. We propose 
to construct a special 50 m radio telescope. Regular timing measurements of about 10 
strong millisecond pulsars will perhaps allow the detection of gravitational waves within 
VLF or at least will give a more stringent upper limits. 

1. Introduction 

Because a background of relic particles gives a snapshot of the state of the uni- 
verse at the time when these particles decoupled from the primordial plasma, relic 
gravitational waves (gravitons) are a potential source of informations on the state 
of the very early universe and the physics at correspondingly high energies, which 
can not be accessed experimentally in any other way 1 . The photons of the cosmic 
microwave background (CMB) give us a snapshot of the state of the universe at 
t ~ 3 x 10 5 years after the big-bang, while the gravitons of the stochastic gravi- 
tational wave background (SGWB) encode in its spectrum the information of the 
universe at t ~ 10 -44 seconds. These properties make the detection of SGWB very 
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interesting, which is the main purpose of the proposal discussed here. 

A SGWB can be characterized by its energy distribution in frequency especially 
the dimcnsionless quantity 

Pcritical 

din/ ' 

where p C riticai = 3c 2 Hq/8ttG is the critical energy density required (today) to close 
the universe (the Hubble constant H = 100/iioo km sec -1 Mpc -1 = 3.2 x 10~ 18 /iioo 
Hz ). Up to now, there have been three observational constraints on SGWB 
spectrum 2 ' 3 : (i) The strongest observational constraint on f2 gw (/) comes from the 
high degree of isotropy observed in the CMB, f2 gw (/) hf 00 < 7xHr n (#o//) for H < 
f < 30H . (ii) The second comes from the standard model of big-bang nucleosyn- 
thesis, ff >10 -8 Hz d\nf f2 gw (/) h\ m < 1CP 5 . (iii) The last observational constraint 
comes from monitoring the radio pulses emitted by a number of stable millisecond 
pulsars for a decade, f2 gw (/ = 10~ 8 Hz) h\ m < 10~ 8 . 

In order to get a lower constraint on SGWB spectrum or even directly detect 
SGWB within the very low frequency band (VLF), 10~ 7 to 10~ 9 Hz, we propose to 
construct a special 50 m radio telescope at National Astronomical Observatories, 
Chinese Academy of Sciences (NAOC) to monitor about 10 strong millisecond pul- 
sars. After reviewing the basics of pulsars as probes for gravitational wave (section 
2), we present our research plan in details (section 3). 

2. Pulsar timing array as efficient gravitational wave detectors 

Pulsars have been believed to be extremely stable natural clocks. Detweiler 4 had 
shown that measurements of pulse arrival time (TOA) may be used to search for 
gravitational waves with a period on the order of years. For SGWB, his found 

243^ 3 / 4 (R*(t)) 
Pgw ~ 208 G • () 

where the residual R(t) is the difference between observed and expected TOA. The 
span of measurements T should be longer than 1/f (with 1// = a few years, the 
measurements may be taken at intervals of say, 10-20 days). 

The observed i?(i) bs consists of deviations caused by various effects. Some of 
them arc systematic, such as the uncertainties in the position and proper motion 
of pulsars, errors in the positions of barycenter. While others are random, for 
example, propagation delays caused by interstellar turbulence, relativistic effects 
by intervening massive objects, irregularities occurred in pulsar parameters, or in 
reference clock. Some of these effects could be modeled or be estimated, and hence 
be corrected. The resulted (R 2 (t)) is therefore the mixture of three components: 
(1) the remaining effects left over after the correction; (2) pulse timing noise which 
is related with (a) signal-to-noise ratio of pulsar observation, (b) pulse shape and 
the definition of arrival time, and (c) method in obtaining the arrival time; (3) 
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background gravitational waves. Thus the quantity (R 2 (t)) from pulsar timing 
observations produces an upper limit on the energy density of SGWB. 

Following Detweilcr's work, a number of observations had been made over a 
decade due to the advent of the discovery of millisecond pulsars. Millisecond pulsars 
have high stability and sharp pulses, therefore the precision of timing measurements 
is extraordinary high, reaching 0.1 /us sometimes. A millisecond pulsar, besides its 
extraordinary high inertia which warrants high stability, the shape of its pulses being 
so sharp has caused a great increase in the precision of timing measurements. As 
a result, the energy density of SGWB in the frequency range around 0.5 cycle/year 
was found to be p < 10~ 6 /? cr iti C ai (see eg. Ref 5). 

This result is encouraging. In particular, a much better sensitivity will be ob- 
tained if one can use two or more pulsars as coincident detectors 6,7 ' 8,9 . Such a 
"pulsar timing array" could be a very efficient SGWB detector. To quantify this, 
consider the Doppler shift from the ith pulsar in the array 9 



where s(t) is the GW signal, on a geometrical factor depending on the line-of-sight 
direction, and rn{t) contains the noise intrinsic to the pulsar timing. The cross- 
correlation between two pulsars then gives 1 



and, since the timing noises from different pulsars are uncorrelated, increasing the 
observation time will limit the cross-correlation towards ctiaj(s 2 {t)). 

Unfortunately, the galactic millisecond pulsars known around 1990, for which 
timing data have been taken already, were concentrated in almost the same region 
of the sky, and correlation analysis of these data does not improve much the result. 
Now the known millisecond pulsars are much more uniformly distributed on the 
sky. It is the right time to construct such a "pulsar timing array" 9 . The Berkeley 
pulsar group led by Don Backer has been working on this direction for a few years. 

3. A special 50 m radio telescope at NAOC 

We propose to construct a special 50m radio telescope at NAOC, devoted to pulsar 
timing measurement. It will work on low frequencies, with 1400 MHz as the highest 
(so that the instruments will not be very expensive) . At least 10 strong millisecond 
pulsars can be chosen from the 27 candidates listed in Table 1 for regular timing 
measurements. Each object should be observed 1 hour every day. With the best 
receiving system today, accumulation of 15 day's measurements will have a sensi- 
tivity of about 0.05 mJy, so that a signal-to-noise ratio above 10 is expected for 
these pulsars. 



Ai/j(t) 



= ais(t) + m(t) 



(3) 



aiaj(s 2 (t)) + ai(sTii) + a j (sn j ) + (r^n,) 



(4) 
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Table 1. 27 candidates to be observed with 50 m telescope (S400 > lOmJy) 



iUOONamc 


Gal-L 


Gal-B 


Period (ms) 


Distancc(kpc) 


c a r\r\ / t \ 

s400(mjy) 


J0034-0534 


111.49 


-68.07 


1.8771818543796 


0.98 


17 


J0613-0200 


210.41 


-9.30 


3.06184403674401 


2.19 


2e+01 


J0751+1807* 


202.73 


+21.09 


3.478770781510 


2.02 


le+01 


J1012+5307 


160.35 


+50.86 


5.25574901198 


0.517 


3c+01 


J1022+1001* 


231.79 


+51.10 


16.452929681440 


0.599 


2e+01 


J1300+1240* 


311.30 


+75.41 


6.2185319388187 


0.624 


2e+01 


J1455-3330 


330.72 


+22.56 


7.987204795504 


0.738 


9.0 


J1623-2631 


350.98 


+15.96 


11.075750876440 


1.8 


2e+01 


J1643-1224 


5.67 


+21.22 


4.62164144656360 


4.86 


8e+01 


J1713+0747* 


28.75 


+25.22 


4.57013652273380 


0.892 


4e+01 


J1730-2304 


3.14 


+6.02 


8.1227979128499 


0.506 


4e+01 


T1 7/1/1 11^/1 
Ji 1 44-1 104 


1 A 70 
14. / y 


+y.io 


4.u / 40400 / oizoyo 


u. loo 


1 Q 
lo 


J1748-2446A 


3.84 


+1.70 


11.563148389585999401 


7.1 




J1804-2717 


3.51 


-2.74 


9.343030681150 


1.17 


15 


J1823-3021A 


2.79 


-7.91 


5.440002278840 


8 


16 


J1824-2452 


7.80 


-5.58 


3.0543146293258 


5.5 


4c+01 


J1857+0943* 


42.29 


+3.06 


5.36210045404154 


1 


3e+01 


J1911-1114 


25.14 


-9.58 


3.6257455713977 


1.59 


31 


J1939+2134* 


57.51 


-0.29 


1.5578064924327 


9.65 


2e+02 


J1955+2908 


65.84 


+0.44 


6.133166488729 


5.39 


2c+01 


J1959+2048 


59.20 


-4.70 


1.60740168480632 


1.53 


2c+01 


J2019+2425* 


64.75 


-6.62 


3.9345240796636 


0.912 


2c+01 


J2051-0827 


39.19 


-30.41 


4.50864174335 


1.28 


2c+01 


J2124-3358 


10.93 


-45.44 


4.93111485914810 


0.248 


17 


J2145-0750 


47.78 


-42.08 


16.0524236584091 


0.5 


100 


J2229+2643* 


87.69 


-26.28 


2.9778192947192 


1.43 


le+01 


J2317+1439* 


91.36 


-42.36 


3.4452510710225 


1.89 


2e+01 



*-Being observed by USA groups 



In addition to what mentioned above, multi-object observation strategy is of 
importance for the detection of SGWB in the following senses: Consider a number of 
pulsars separated at different parts in the sky, gravitational wave from any direction 
sweeping across the Earth will be marked by the change of the rates of local clock at a 
frequency /. Owing to the quardrupole character of the gravitational radiation, such 
effect will exactly in two pulsar data if their directions are 180° apart. So that R(t), 
in the form of random fluctuations, contained in the record of timing measurements 
of a pulsar would have an exactly similar shape with the R(t) record of a pulsar. As 
the period of the gravitational wave 1// being very long, observations of two pulsars 
within a time interval of one day can be thought as nearly simultaneous. R(t) \> s 
of the two records, when cross-correlated, will reveal the existence quantitatively 
of the background gravitational waves. In practice, we will have a number of pairs 
of pulsars at different separation angles, by the same principle, cross-correlations 
performed to the i?(i) bs of all the pairs will produce a better result than that from 
a single pair. 

Finally, we would like to remark that, among these sources of perturbations 
contained in R(t), many of them are of interest in their own right. The timing 
measurement data of these objects may be used to, for examples, (l)improve the 



Special Purpose Pulsar Telescope 



cphcmcris, (2)monitor the intrinsic perturbations of pulsars, (3)establish a time 
standard for long time scales. 
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